Percolation Critical Exponents in Scale-Free Networks 



Reuven Cohen 1 0, Daniel ben-Avraham 2 , and Shlomo Havlin 1 

1 Minerva Center and Department of Physics, Bar-Ban university, Ramat-Gan, Israel 
2 Department of Physics, Clarkson University, Potsdam NY 13699-5820, USA 



(N 
O 
O 
(N 

00 



We study the behavior of scale-free networks, having connectivity distribution P(k) ~ k~ x , close 
to the percolation threshold. We show that for networks with 3 < A < 4, known to undergo a 
transition at a finite threshold of dilution, the critical exponents are different than the expected 
mean-field values of regular percolation in infinite dimensions. Networks with 2 < A < 3 possess 
only a percolative phase. Nevertheless, we show that in this case percolation critical exponents are 
well defined, near the limit of extreme dilution (where all sites are removed), and that also then the 
exponents bear a strong A-dependence. The regular mean-field values are recovered only for A > 4. 
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Large networks have been attracting considerable in- 
terest |],pL[|[|j|,@] . A useful characterization of such net- 
works is their degree distribution, P(k), or the probabil- 
ity that an arbitrary node be connected to exactly k other 
nodes. Many naturally occurring networks (the Internet, 
nets of social contacts, scientific collaborations, ecologi- 
cal nets of predator-pray, etc.,) exhibit a power-law, or 
scale-free degree distribution: 



P{k) = ck~ 



k > m 



(1) 



where m is the minimal connectivity (usually taken to be 
m = 1) and c is a normalization factor. 

Imagine that a large network is diluted, by random 
removal of a fraction p of its nodes. If p is small, one 
expects that the network remains essentially connected: 
there exists a giant component of connected nodes that 
constitutes a finite fraction of the total size of the original 
network. Above a certain threshold of dilution, p c , the 
giant component disappears and the network effectively 
disintegrates. This is the problem of percolation |7],^] as 
defined on networks. (Percolation was originally defined 
on regular lattices.) The giant component corresponds 
to the infinite incipient cluster or spanning cluster which 
forms only in the percolating phase, at dilutions smaller 
than the percolation threshold p c . 

Percolation in scale-free networks is widely recognized 
as a key problem of interest [[l]]|,||,||,|],D . Applications 
range from the robustness of communication networks 
(the Internet, phone networks) in the face of random fail- 
ure (removal) of a fraction of their nodes (e.g., random 
breakdown of routers), to the efficiency of inoculation 
strategies against the spread of disease, to the ecological 
impact of the extinction of some species. 

For a random network of arbitrary degree distribu- 
tion, the condition for the existence of a spanning cluster 
is SMajIl: 



(k) 



(2) 



Suppose a fraction p of the nodes (and their links) are 
removed from the network. (Alternatively, a fraction 
q = 1— p of the nodes are retained.) The original degree 
distribution, P(k), becomes 



CXJ j 

p'(fc) = £p(0(J(i 



(3) 



l=k 



In view of this new P'(k), Eq. (|2|) yields the percolation 
threshold: 



<7c = 1 - Pc 



1 



K — 1 



(4) 



where k is computed with respect to the original distri- 
bution, P(k), before dilution. 

Applying the criterion (^) to scale-free networks one 
concludes that for A > 3 a phase transition exists at a 
finite q c , whereas for 2 < A < 3 the transition takes place 
only at the extreme limit of dilution of q c = |||| . Here 
we concern ourselves with the critical exponents associ- 
ated with the percolation transition in scale-free networs 
of A > 2. We show that critical exponents are well de- 
fined (and we compute their A-dependent values) even in 
the non-percolative regime of 2 < A < 3. In the per- 
colative range of A > 3, we find A-dependent exponents 
which differ from the regular mean-field values of per- 
colation in infinite dimensions. This irregularity can be 
traced to the fat tail of the distribution (1). The regular 
mean-field exponents are recovered only for A > 4. 

In UJl^] a generating function is built for the connec- 
tivity distribution: 
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G (x) = J2 p ( k > k - 



(5) 



fe=0 



The probability of reaching a site with connectivity k by 
following a specific link is kP(k)/(k) p|]9l,|l2|| , and the 
corresponding generating function for those probabilities 
is 



£fcP(fc); 



Jfc-1 



EkP(k) 



dx 



G (x)/(k) 



(6) 



Let Hi(x) be the generating function for the probability 
of reaching a branch of a given size by following a link. 
After a dilution of a fraction p of the sites, Hi (x) satisfies 
the self-consistent equation 



H^x) = l-q + qxG x (Hi{x)) 



(7) 



Since Gq(x) is the generating function for the connectiv- 
ity of a site, the generating function for the probability 
of a site to belong to an n-site cluster is 



H (x) = l-q + qxG (H 1 (x)) 



(8) 



-flo(l) is the probability that a site belongs to a cluster 
of any finite size. Thus, below the percolation transition 
Hq(1) = 1, while above the transition there is a finite 
probability that a site belongs to the infinite spanning 
cluster: P^ = 1 — Hq(1). It follows that 



P 00 (q) = q(l-J2P(k)u k ) 



(9) 



k=0 



where u = -Hi(l) is the smallest positive root of 



.k-l 



(10) 



This equation can be solved numerically and the solu- 
tion may be substituted into Eq. (|9|), yielding the size 
of the spanning cluster in a network of arbitrary degree 
distribution, at dilution q 

We now compute the order parameter critical expo- 
nent j3. Near criticality the probability of belonging to 
the spanning cluster behaves as ~ (q — q c )@ . For 
infinite-dimensional systems (such as a Cayley tree) it is 
known that /3 = 1 [0J|Jl^]. This regular mean- field re- 
sult is not always valid, however, for scale-free networks. 
Eq. (||) has no special behavior at q — q c ; the singular 
behavior comes from u. Also, at criticality Poo = and 
Eq. (JsJ) imply that u = 1. We therefore examine Eq. jlO| ) 
for u = 1 — e and q = q c + 5: 



l-e = l-q c 



(gc + S) 
(k) 



oo 



k=0 



kP{k){l~e) k - 1 . (11) 



j^feP(Ar)«* -1 ~ (k) - {k(k - l))e 



k=0 

+^(k(k -l)(k- 2))e 2 + ■ • • + cr(2 - A)e A - 2 , (12) 

where the highest-order analytic term is 0(e n ), n = 
[A — 2J . Using this in Eq. Q, with q c = 1/(k - 1) = 
(k)/(k(k-l)), we get 

iHk ~ ) 1))2 S= 1 -(k(k-l)(k-2))e + ... 



-cT{2 - A)e 



A-3 



(13) 



The divergence of 8 as A < 3 confirms the lack of a phase 
transition in that regime. Thus, limiting ourselves to 
A > 3, and keeping only the dominant term as e — > 0, 
Eq. (O) implies 




c(fe)r(2-A) 
2(fc(fc-l)) 2 



(fe)(fe(fe-l)(fe-2)) ' 



5— 3<A<4, 



A > 4. 



(14) 



Returning to Poo, Eq. (||), we see that the singular con- 
tribution in e is dominant only for the irrelevant range 
of A < 2. For A > 3, we find P^o — q c {k)e ~ (q - q c )$ . 
Comparing this to ( |l4| ) we finally obtain 







5^ 3<A<4, 
1 A > 4. 



(15) 



We see that the order parameter exponent /3 attains 
its usual mean-field value only for A > 4. Moreover, for 
A < 4 the percolation transition is higher than 2nd-order: 



for 3 + zr-j < A < 3 



ig the transition is of the nth- 



order. The result (15) has been reported before in [||, 
and also found independently in a different but related 
model of virus spreading |lj,|l^|. The existence of an 
infinite-order phase transition at A = 3 for growing net- 
works of the Albert-Barabasi model, has been reported 
elsewhere |l6|,[l7j ■ These examples suggest that the criti- 
cal exponents are not model-dependent but depend only 
on A. 

For networks with A < 3 the transition still exists, 
though at a vanishing threshold, q c = 0. The sum in 
Eq. ( fll| ) becomes: 



kP(k)u k - 1 - (k) + cT(2 - A)e 



A-2 



(16) 



k=0 



Using this in conjunction with Eq. (|To|), and remember- 
ing that here q c = and therefore q — S, leads to 



'-er(2-A)\5±x._i 
' = 1 (k) ) S3 - X > 



(17) 



The sum in (Jill) has the asymptotic form 



which implies 



2 



A 



2 < A < 3 . 



(18) 



In other words, the transition in 2 < A < 3 is a mirror 
image of the transition in 3 < A < 4. An important dif- 
ference is that q c = is not A-dependent in 2 < A < 3, 
and the amplitude of Poo diverges as A — > 2 (but remains 
finite as A — > 4). 

In it was shown that for a random graph of ar- 
bitrary degree distribution the finite clusters follow the 
usual scaling form: 



s e ' 



(19) 



Here s is the cluster size, and n s is the number of clus- 
ters of size s. At criticality s* ~ \q — q c \~ a diverges and 
the tail of the distribution behaves as a power law. We 
now derive the exponent r. The probability that a site 



belongs to an s-cluster is p s = sn s 
ated by H : 



PsX 



and is gener- 



(20) 



The singular behavior of Hq(x) stems from H\(x), as can 
be seen from Eq. (^) . Hi (x) itself can be expanded from 
Eq. (0), by using the asymptotic form ( |l2|) of G\. We 
let x = 1 — e, as before, but work at the critical point, 
q = q c . With the notation </>(e) = 1 — Hi(l — e), wc finally 
get (note that at criticality -Hi(l) = 1): 



(k(k 



(1- 
l)(fc 



9c 



2)} 



2(k) 



+ he 



r(2-A) 

(k) ' 



,A-2i 



(21) 



From this relation we extract the singular behavior of Hq: 
4> ~ e y . Then, using Tauberian theorems |Q it follows 
that p s ~ s _1_y , hence r = 2 + y. 

For A > 4 the term proportional to 4> x ~ 2 in ( pi] ) may 
be neglected. The linear term e</> may be neglected as 
well, due to the factor e. This leads to <j) ~ e 1 / 2 and to 
the usual mean-field result 



5 

T = — 

2 



A > 4 



For A < 4, the terms proportional to ec; 
glected, leading to 4> ~ e 1 /( A ~ 2 ) and 



(22) 



may be ne- 



1 



A 



2A-3 
A-2 



2 < A < 4 



(23) 



Note that for 2 < A < 3 the percolation threshold is 
strictly q c = 0. In that case we work at q = 5 small 
but fixed, taking the limit S — * at the very end p9| . 
For growing networks of the Albert-Barabasi model with 
A = 3, it has been shown that sn s cx (slns)~ 2 ^7^. This 
is consistent with r = 3 plus a logarithmic correction. 
Related results for scale free trees have been presented 
in J2l. 



The critical exponent a, for the cutoff cluster size, may 
be too derived directly. Finite-size scaling arguments pre- 
dict that 



q c {oo)-q c {N) - N- 



= N~ 



(24) 



where N is the number of sites in the network, v is the 
correlation length critical exponent: £ ~ (q — q c ) v , and 
d is the dimensionality of the embedding space. Using 
a continuous approximation of the distribution (]l|) one 
obtains lb| 



2-AA K 3 - x -m 3 - x 



A J K 2 -> 



,2-A 



(25) 



where K ~ jSf 1 /^ : ) is the largest site connectivity of 
the network. For 3 < A < 4, this and Eq. (0) yield 



q c (oo)~q c (N) - Ak ~ K 



3-A 



N 1 



which in conjunction with Eq. ( 
A-3 



A-2 



leads to 
3 < A < 4 . 



(26) 



(27) 



For A > 4 we recover the regular mean-field result 
a = 1/2. Note that Eqs. (fj), (pf), © are consistent 
with the known scaling relation: ct/3 = r — 2 ]7],p|,Jl3[ . 



For 2 < A < 3, g c (oo) 
iV (A-3)/(A-i) > Therefore 

3-A 



A-2 



and q c (N) 



2 < A < 3 , 



(28) 
2 (cf 



again consistent with the scaling relation aj3 = r 
Eq. (g)). 

Through similar scaling relations, any two of the perco- 
lation exponents discussed above determine the remain- 
ing (static) percolation exponents. Thus, for example, 
the exponent 7, governing the average size of finite clus- 
- <7c) 7 , follows from the scaling relation 
In view of Eqs. @, @, (M), we get 



ters: (s) ~ (q 
7 = (3 - r)ja 



7 



A > 3, 
2 < A < 3. 



(29) 



The negative value of 7 in the range 2 < A < 3 is, 
at first sight, surprising. However, recall that for that 
range only the percolating phase exists. At the transi- 
tion point, q c — 0, all clusters have size zero. As q in- 
creases, finite-size clusters begin to show, consistent with 
(a) ~ (q- q c )~ 7 - q. 

In summary, we have shown that percolation critical 
exponents in scale-free networks bear a strong depen- 
dence upon the degree distribution exponent A. This 
is true even in the range 3 < A < 4, where the percola- 
tion transition occurs at a finite threshold q c : the regular 
mean-field behavior of percolation in infinite dimensions 
is recovered only for A > 4. Moreover, critical exponents 



3 



are well defined also in the most physically relevant re- 
gion of 2 < A < 3 (the case of most naturally occurring 
networks), despite the lack of a non-percolating phase. 
In this regime too the critical exponents depend strongly 
on A. 

An interesting exception is the exponent 7, which at- 
tains constant values independent of A, Eq. (^). This 
suggests some underlying deeper principle which we were 
not yet able to uncover. 
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